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Recent experiments show that double layer quantum Hall systems may have a ground state with 
canted antiferromagnetic order. In the experimentally accessible vicinity of a quantum critical point, 
the order vanishes at a temperature Tkt = kH, where H is the magnetic field and k is a universal 
number determined by the interactions and Berry phases of the thermal excitations. We present 
quantum Monte Carlo simulations on a model spin system which support the universality of k and 
determine its numerical value. This allows experimental tests of an intrinsically quantum-mechanical 
universal quantity, which is not also a property of a higher dimensional classical critical point. 



Quantum Hall systems offer attractive, tunable labora- 
tories for investigating zero temperature quantum transi- 
tions between states with different spin magnetizationsa. 
Their ground states are determined almost entirely by 
the Coulomb interactions between the electrons, and the 
typical Coulomb exchange energy is usually much larger 
than the Zeeman energy in the external field; as a result, 
the spins are often not fully polarized in the direction 
of the applied field, and can realize different magnetic 
configurations which optimize the Coulonab interactions. 

Three separate recent experimentsacl have studied 
magnetic transitions in bilayer quantum Hall systems at 
total filling fraction v = 2. When the layers are well sep- 
arated, each layer forms a fully polarized, ferromagnetic 
state with all states in the lowest Landau level occupied. 
The parallel alignment of all spins is induced mainly by 
an intralayer ferromagnetic exchange interactionO, but is 
also compatible with the Zeeman coupling which orients 
the spins in the direction of the magnetic field. When 
the layers are closer to each other, there is a significant 
interlayer antiferromagnetic exchange interactions which 
eventually prefers a spin singlet ground state. The transi- 
tion from the fully polarized ferromagnet to the spin sin- 
glet quantum paramagnet could, in principle, be a direct 
first-order transition; however, it was theoreticallyLl and 
experimentallyBBI found to occur via a softening of the 
energy of a single spin-flip excitation in the ferromagnet, 
suggesting an intermediate phase with canted spin or- 
dering!!], bounded by second-order transitions. Detailed 
theoretical predictions of the phase diagram have been 
madeO'E2l, and are in good agreement with recent light 
scattering experimentsB. 

In this paper, we shall use quantum Monte Carlo sim- 
ulations to study a bilayer quantum spin model which 
has a phase diagram closely related to that of the v = 2 
bilayer quantum Hall system. In particular, all quantum 
and nonzero temperature (T) phase transitions of the 



two models are expected to be in the same universality 
classes. Our focus will be on the phase with canted spin 
ordering. This ordering breaks a U(l) spin rotation sym- 
metry, and it is expected that the symmetry will be re- 
stored at nonzero T by a Kosterlitz-Thouless (KT) phase 
transition at T = Tkt- In general, the value of Tkt de- 
pends upon microscopic details of the Hamiltonian, but 
in the vicinity of a certain quantum— critical point (see 
discussion below and Fig. pi) it obeysO 
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where H is the external magnetic field (the electron gy- 
romagnetic ratio and the Bohr magneton have been ab- 
sorbed into the definition of H) and k is a non-trivial 
universal number. It turns out that our quantum spin 
model has two separate quantum critical points for which 
(1) is expected to be valid. Our simulations verify that 
(I) is indeed obeyed near both critical points; moreover, 
the values of k determined at the two points are identical 
to within the numerical accuracy, and this supports the 
claimed universality of k. The same value of k should also 
apply to the bilayer quantum Hall systems, and this is 
a quantitative theoretical prediction which can be tested 
in future experiments. 

The bilayer quantum spin model has the Hamiltonian 



H = Y J • S 2l - H • (S u + S 2 

i 

<ij> 



(2) 



where S a i are quantum spin- 1/2 operators in 'layers' 
a = 1,2 residing on the sites, i, of a two-dimensional 
square lattice, H = (0, 0, H) is the external magnetic 
field, and J±, J are intra- and inter-layer exchange con- 
stants respectively. We will take J± > antiferromag- 
netic, but allow J to take either sign. 
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FIG. 1. Ground states of TL. The arrows denote the mean 
orientations of the spins in the two layers, for the case where 
H points vertically upwards; in the SS phase the spins have no 
definite orientation. The in-plane ordering wavevectors of the 
x, y components of the spins in the C phase are indicated. The 
boundaries of the SS phase approach the points M M (/j, — 1, 2) 
as H/J± ~ | to — WuV where w = J/J±, Wi = 0.398 (also 
determined in Refs. [l^,[13p, u>2 = —0.435 (also determined in 
Ref. |l4j), and v ~ 0.7 is the correlation length exponent of the 
three dimensional classical Heisenberg ferromagnet. There is 
a Kosterlitz Thouless transition at non-zero T above the C 
phases, and the result ([j|) applies above the shaded regions. 



The., model TL has been studied intensively in recent 
yearal3il3 for the case H = 0. Using the methods of 
Ref. [l5], and numerical and exact analytical results to be 
discussed below, we obtained the H ^ phase diagram 
shown in Fig. ^. We will now discuss phases in Fig. ^ in 
turn, and also indicate the nature of the quantum tran- 
sitions between them: 

(i) Fully Polarized Ferromagnet (FPF): For large enough 
H, the exact ground state is simply the state with all 
spins up i.e. in the eigenstate of S a i Z with eigenvalue 
1/2. The first excited state is a single spin- flip, and its 
excitation energy can be determined exactly: at momen- 
tum k = (k x , k y ) it is = H — J± — J(2 — cos k x — cos k y ). 
For J > ( J < 0) this has a minimum at k = (it, ir) 
(k = (0,0)). Stability of the FPF state requires that 
f£ => 0, and the point where the minimum energy first 
vanishes exactly determines the FPF phase boundary 
shown in Fig. [y. The single spin-flips Bose condense at 
this boundary, leading to a phase with canted spin order 
to be described below. This transition is in the universal^ 
ity class of the dilute Bose gas quantum critical pointed 
with dynamic exponent z — 2. 

(ii) Spin Singlet (SS): This is the spin singlet quantum 
paramagnet with no broken symmetries and a gap to all 
excitations. The exact ground state is known only for 
J = 0, when the system decouples into pairs of spins an- 
tiferromagnetically coupled by Jj_, which therefore form 
a spin singlet valence bond (see Fig. P . The ground state 
for J 7^ is adiabatically connected to this decoupled 



state, and its wavefunction can be determinedlll in an ex- 
pansion in powers of J/ J±. For H = 0, the lowest excited 
state is triplet particle, again adiabatically connected to. 
the J — limit, whose dispersion has been computed^ 
in a series in J/J±. Turning on a nonzero H leads to no 
change in the wavefunctions, but the energy of the par- 



ticle changes by e? 



mH , with in = 1,0, —1 the S z 



quantum number. As for the FPF phase, the boundary 
of the SS phase is the line where the minimum excitation 
energy vanishes. For J > 0, the boundary is at H/ J± = 
l-2w- (3/2)w 3 - (3/2)w 4 -I- 0(w 5 ), with w = J/J±, 
while for J < it is at H/Jj_ = l + 2w-{7/A)w 4 + 0{w 5 ). 
For H > 0, there is a Bose condensation of the particle 
at this line, which is in the same universality class as the 
boundary of the FPF phase, and also leads to the same 
canted phase. Precisely at H = (the points Mi, Mi in 
Fig. [l]) the nature of the transition is different, and will 
be discussed below. 

(Hi) Canted (C): TL has a symmetry of rotations about 
the axis (z) of the applied field, and this is broken in this 
phase. The spin operators have the expectation values 
(S lz ) = (S 2z ) and (Si x , y ) = -{S 2x , y ) + 0. The 
z expectation values are independent of i, while the x,y 
expectation values have staggered (uniform) arrangement 
on the two sublattices with each layer for J > (J < 0). 
Associated with the broken symmetry, there is linearly 
dispersing Goldstone mode of spin-wave excitations cor- 
responding to slow rotations of the order parameter in 
the x, y plane. 

(iv) Neel (N): This is reached in the H = limit of the C 
phase, when (S az ) — 0. Now H has the full SU(2) spin 
rotation symmetry, and so the spin expectation values in 
the x, y plane can actually point along any direction in 
spin space. 

The vicinities of the critical points M^2 are of par- 
ticular,-] interest in this paper. Here the system is ex- 
pectecl3't2l to be described by a continuum quantum 
field theory which can be expressed in terms of a unit 
length field n(x, r) (r is imaginary time). For Mi, 
n « E^M^+MSi* - S w ) and for M 2 , n cx 

SieAkC^it - S 2 ;), where N x is an averaging neighbor- 
hood of x. The field theory has the action 
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with c a velocity and g a coupling constant which tunes 
the value of J/J±. At H = this theory can be 
reinterpreted as the real partition function of a three- 
dimensional classical Heisenberg ferromagnet at non-zero 
temperature. However no such classical interpretation is 
possible for H =/= 0: notice then that the action is com- 
plex (even in imaginary time), as it includes the Berry 
phase of the precession of the n quanta about the ap- 
plied field. Applying a field at the scale-invariant critical 
points M\ t 2, we can conclude that all characteristic tepr 
peratures will be determined by H . Scaling argumentsliil, 



2 



which rely on the fact that H appears in S as the time 
component of a 0(3) non-Abelian gauge field (which is in 
turn related to the fact that H couples to the conserved 
total spin), show that H scales as an inverse time; as tem- 
peratures also scales as inverse time, we are then lead to 
(|l|) for the critical temperature at which the field-induced 
canted order will disappear. 

We turn now to our quantum Monte CapLo results, ob- 
tained using the powerful loop algorithmic successfully 
used in-cecent studies of quantum HeisenbergEZI and XY 
modelsEa. All weights are positive in this basis if we 
apply the field along the axis of quantization, and the 
Berry phases have been transformcd-4nto the quantiza- 
tion of the spins on the world lineali. The loop algo- 
rithm slows down severely for H ^ Olj, as a loop that 
changes the magnetization by AM picks up a flipping 
probability exp(HAM/T); this leads to an exponential 
increase of the autocorrelation times from T lnt < 1 to 
Ti n t °c exp(H/T). Fortunately it is sufficient to simulate 
at not too low T, where H/T < 4 and r int < 100, but 
our systepi sizes, L, are not as large as in earlier H = 
studiesOlia. 

We obtained Tkt following the method of Harada and 
Kawashimalia for the quantum XY model. The spin stiff- 
ness, p s , of the [7(1) ordering in the x — y plane was ob- 
tained from p s = T(Wj + Wy) /2, where W x , y are the 
total winding numbers in the two directions. The im- 
proved estimators in Ref. [l8] have to be modified, since 
the flipping probabilities of loops are no longer all equal, 
and a multi-loop algorithm is necessary. For L = oo, 
p s is finite for T < Tkt, continuously decreases to the 
universal value p s {T — Tkt) — {^/^)Tkt at Tkt and is 
zero for all T > Tkt- However, for L finite, p s is nonzero 
at all T, and is expected to converge to the Lu^-po limit 
with the finite size scaling form at T — Trt^'E^- 

n Ps /2T = 1 + [2 log(L/L (T))]-\ (4) 

Hence, good estimates for Tkt can be obtained by plot- 
ting l/(irp s /T — 2) — logi as a function of L. As 
L is increased this quantity converges to the constant 
(— logLo) at Tkt, and diverges to ±oo for T > Tkt and 
T < Tkt respectively. Our data (a representative exam- 
ple is shown in Fig. ||) are clearly consistent with these 
expectations. Note that we obtained-Xo ~ 5 to 10, com- 
pared to La = 0.23 in the XY modeltS; so we cannot use 
the more elaborate fitting techniques used in Ref. [l^, as 
we are not deep enough in the asymptotic scaling regime. 

First, we checked for a KT transition for the decou- 
pled single-layer case (Jj_ = 0). Our results for p s 
(per layer) are shown in Fig. [3|a. We find KT tran- 
sitions at T KT I J = 0.190(5) for H/J = 0.1 and at 
Tkt/ J = 0.215(5) for H/J = 0.2. Surprisingly however, 
anomalous non-monotonic finite size scaling behavior was 
found below Tkt- As can be seen in Fig. ||a, when in- 
creasing L, first p s decreases, then increases again, and 
is asymptotically expected to decrease again, like in the 
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FIG. 2. Scaling plot of l/(7rp s /T — 2) — log L versus L. This 
quantity, which diverges to +oo for T > Tkt and to — oo for 
T < Tkt allows a reliable estimate of Tkt- 
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FIG. 3. Spin stiffness p s in a magnetic field H — 0.2J of 
TL for (a) the single layer model (J± = 0) and (b) vertically 
above the critical point Mi of Fig. [j], as a function of T for 
various L. The dashed lines are the lines of universal values 
at T KT : (tt/2) Ps (T kt ) = Tkt- 
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the universal temperature of a KT transition in the vicin- 
ity of a quantum critical point. This universal quantity 
relies on an underlying interacting quantum field theory 
with complex Berry phases in d = 2, and experiments 
to measure it in bilayer quantum Hall systcmjL.will be of 
considerable interest. Current measurements™ of Tkt 
are in the vicinity of ([!]) , but more detailed measurements 
in the shaded region of Fig [l] axe necessary, possibly by 
pressure tuning of the g-factorcil 
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University of Tokyo, using a parallelizing Monte Carlo li- 
brary in C++ developed by one usEJ. 
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FIG. 4. The critical temperature Tkt as a function of field 
H for the square lattice (isolated layers) and for the antifer- 
romagnetic (AFM) and ferromagnetic (FM) bilayer models 
close to both critical points. 



XY model. A similar anomalous scaling was found in 
this model by Lavalle et aZo for the ground state energy 
in subspaces of nonzero spin S oc L? . 

We now turn to fields above the critical points M\ and 
M2. We performed simulations in fields H/J = 0.1, 0.15, 
0.2, 0.4, 0.6 and 1. Here no anomalous finite size scaling 
is observed, as can be seen in Fig. |3|b (the p s now is 
the total value, not per layer). We can thus confidently 
use the finite size scaling analysis to determine Tkt as a 
function of H, and show the results in Fig. 

We fit the results of Fig. § H/\ J\ < 0.2 to T KT = 
a(6) + kH, where the offset a(6) accounts for the error in 
our determination of the positions of Mi )2 . We obtained 
for the slope k = 0.35(7) at M x and k = 0.37(5) at M 2 . 
These numbers agree very well. The error bars are our 
estimate of strict upper and lower bounds, and not one- 
sigma confidence intervals. We obtained a(6) « 0.01J, 
comparable with what the H — 0, T = p s could be, 
given the uncertainty of about 0.4% in the determina- 
tion of the critical points (compare the results for slightly 
different coupling ratios in Fig. |]). 

Corrections to scaling from not being in the contin- 
uum limit appear at H ~ J and lead to Tkt — kH(1 — 
bH/\J\): these are smaller and were ignored in above fits 
for H/\J\ < 0.2. Using above estimates for k and a(S) 
and a small correction b « 0.3 we can however fit all 
our results up to H = 0.6|J|. This correction slightly 
increases our final combined estimate: 



kw 0.38 ±0.06; (5) 

we aJsq recall t he leading result in the e = 3 — d expan- 
sioifM k = ^33/107r 2 e w 0.58. 

To conclude, we have obtained a numerical estimate for 



Perspectives in Quantum Hall Effects, S. Das Sarma and 
A. Pinczuk eds, Wiley, New York (1997). 

2 V. Pellegrini et al Science 281, Aug 7 (1998). 

3 J.G.S. Lok et al, |sond-mat/9804256| . 

4 A. Sawada et al, Phys. Rev. Lett. 80, 4534 (1998). 

5 Y.A. Bychkov, S.V. Iordanskii and G.M. Eliashberg, JETP 
Lett. 33, 143 (1981); C. Kallin and B.I. Halperin, Phys. 
Rev. B 30, 5655 (1984). 

6 S. Das Sarma, S. Sachdev and L. Zheng, Phys. Rev. Lett. 
79, 917 (1997). 

7 L. Zheng, R.J. Radtke and S. Das Sarma, Phys. Rev. Lett. 
78, 2453 (1997). 

8 A. Pinczuk et al. Bull. Am. Phys. Soc. 41, 482 (1996); A.S. 
Plaut, ibid. 41, 590 (1996). 

9 V. Pellegrini et al. Phys. Rev. Lett. 78, 310 (1997). 

10 S. Das Sarma, S. Sachdev and L. Zhe ng, Phys. Rev. B 58, 
August 15 (1998); |cond-mat/97093l4 

11 S. Sachdev, Z. Phys. B 94, 469 (1994). 

12 T. Matsuda and K. Hida, J. Phys. Soc. Jpn 59, 2223 
(1990); K. Hida, ibid 59, 2230 (1990); A.J. Millis and 
H. Monien, Phys. Rev. Lett. 70, 2810 (1993); A.W. Sand- 
vik and D.J. Scalapino, ibid 72, 2777 (1994); A.W. Sand- 
vik, A.V. Chubukov and S. Sachdev, Phys. Rev. B 51, 
16483 (1995); C.N.A. van Duin and J. Zaanen, Phys. Rev . 



Lett. 78, 3019 (1997); Zheng Weihong, |eond-mat/9701214| ; 
M.P. Gelfand, Zheng Weihong, C.J. Hamer and J. Oitmaa, 
Phys. Rev. B 57, 392 (1998); V.N. Kotov, O. Sushkov, 
Zheng Weihong and J. Oitmaa, Phy s. Rev. Lett. 80, 5790 
(1998); K. Hida, fcond-mat/9712293|; L. Yin, M. Troyer and 



S. Chakravarty, Europhys. Lett. 44, 559 (1998). 

13 M.P. Gelfand, Phys. Rev. B 53, 11309 (1996). 

14 Y. Matsushita et al, J. Phys. Soc. Jpn 66, 3648 (1997). 

15 S. Sachdev and T. Senthil, Annals of Phys. 251, 76 (1996). 

16 H. G. Evertz et al. Phys. Rev. Lett. 70, 875 (1993); B. B. 
Beard and U.-J. Wiese, Phys. Rev. Lett. 77, 5130 (1996). 

17 M. Troyer et al, Phys. Rev. Lett. 76, 3822 (1996); M. 
Troyer et al, Phys. Rev. B 55, R6117 (1997); M. Troyer et 
al, J. Phys. Soc. Jpn. 66, 2957 (1997). 



4 



K. Harada and N 
(1997) 



Kawash ima, Phys. Rev. B 55, R11949 



ond-mat/980309C. 



' For H 7^ 0, any basis will have either complex Berry phases 
or quantized spins. In contrast, for H — 0, 5 can be repre- 
sented by a continuous sp in classical ferromag net in d = 3. 

1 V. A. Kashurnikov et of., |cond-mat/9802294 



H. Weber and P. Minnhagen, Phys. R ev. B 37, 5986(198 7). 

C. Lavalle, S. Sorell a and A. Parola, :gnd-mat/9709174 . 

D. R. Leadley et al, 3ond-mat/9805357. 

M. Troyer, B. Amnion and E. Heeb, Lecture Notes in Com- 
puter Science (Springer Verlag, in press). 



5 



